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Abstract 

We study the doubly charmed baryonic B decays B — > S C A C and B — ► A^A~K, recently observed 
by BELLE. We find that the unexpected large branching ratios (BRs) could be ascribed to the final 
state interactions (FSIs), which are dictated by B -> DD+ — >■ H C A+ and 5 — ► D°D°K — >■ A+A~if. 
By utilizing the same mechanism, we predict that the BRs for 5+ E°£+ and -> Y,fk±K + 
decays could be as large as BR{B + — > 3°A+) and 10 -4 , respectively. In addition, extending the 
FSIs to the processes associated with the creation of the ss pair, BR(B + — ► ^ c £<^) at percent level 
is achievable. 
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There have been many baryonic B decays studied at B factories. They are the charmless 
decays B — ► p(A)p(A)K , B — > pAn jlj and B — > p(A)p(A) in which the three-body 
baryonic decays have branching ratios (BRs) in mag nitude of O(10~ 6 ) while the BRs of two- 
body decays are limited to be less than 6.9 x 10~ 7 . In addition, the single charmed baryonic 
decay B° -> ^~pir + is measured to be (2.38+°;| ± 0.41 ± 0.62) x 10~ 4 ^ while B° -> A+p 
and 5" -> S c (2455)°p are (2. 19^ ±0.32 ±0.57) x 10~ 5 and (3.67^ ±0.36 ±0.95) x 10~ 5 
jj], respectively. Phenomenogically, the BRs of three-body baryonic B decays are roughly 
one order of magnitude larger than those of corresponding two-body decays. The preference 
could be understood by the threshold enhancements occurring in the near threshold of 
baryon-pair invariant mass QQ. Moreover, one can easily and that in both three-body and 
two-body decays, the BRs for single charmed baryonic processes are two orders of magnitude 
larger than those for charmless decays. The differences could be attributed to the Cabibbo- 
Kobayashi-Maskawa (CKM) []| matrix elements, where single charmed (charmless) processes 
are associated with V C b(V u b). A more detailed review could be referred to Ref. 8|. 

However, when BELLE observes the two charmed baryons in the final state, the char- 
acters appearing in the single charmed and charmless decays subsequently are changed. 
According to BELLE's results, the BRs BR(B + -> A+A~K + ) = (6.5iJ;§ ± 1.1 ±3.4) x 10~ 4 
and BR(B° -> A+A~K°) = (7.9±|;f ± 1.2 ± 4.1) x 10~ 4 are measured with statistical sig- 
nificance of 16.4cr and 6.6cr, and the products of BRs for two-body decays are BR(B + — > 
2°A+) x BR(E° C -> S+tT) = (4.8^;° ± 1.1 ± 1.2) x 10~ 5 and BR(B° -> E~A+) x BR(E~ -> 
S + 7T _ 7T _ ) = (9.3i|g ± 1.9 ± 2.4) x 10^ 5 with 8.7a and 3.8a significance, respectively. If we 
use the theoretical calculations BR(E° C -> E~ir + ) = 1.3% and BR(E+ -> S°vr + ) = 3.9% 
11| and the data BR(E+ -> E°ir + )/BR(E+ -> E-tt + tt + ) = 0.55 ± 0.16 Q, we can get 
B + -> H°A+ « 4.8 x 10~ 3 and B° -> H"A+ « 1.2 x 10" 3 Q. Based on these results, one 



can find that there definitely appear two puzzles: (I) in terms of the BR of B° — > T,~~pir + , 
one can immediately find that due to phase space suppression, the BR of B + — > A+A~K + is 
O(10~ 6 ) which is two orders of magnitude smaller than the observation; (II) It is expected 
that BR(B — > E C A+) ~ BR(B° — > A~p) in theoretical calculations [13|, but the reality 
chooses BR(B — > H C A+) >> BR(B° — > A~p). In this paper, we are going to investigate 
the possible mechanism to solve the unexpected large BRs in doubly charmed baryonic B 
decays. 

It has been noticed that final state interactions (FSIs) may play an important role on the 
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BRs, CP asymmetries (CPAs) and polarizations of vector mesons in two-body charmless and 



charmed mesonic B decays 



14j . According to our previous discussions, the inconsistency in 



two charmed baryons production results from ordinary estimations and QCD calculations. 
That is, it should exist a new mechanism to enhance the BRs for the decays B — > S r A+ and 

n 

B —> A+A C K. In Ref. |6|, two different mechanisms have been introduced, in which the 
two-body and three-body decays are governed by the a, 7r° and 7r =l= meson exchanges and 
charmonium-like resonance X c5 , respectively. However, according to the results of BELLE 
jjj, there is no evidence for the existence of charmonium-like resonance. In addition, since 
the doubly charmed baryonic B decays are found in two-body and three-body decays si- 
multaneously, we speculate that they are induced by the similar mechanisms. Hence, we 
consider that there exist other mechanisms which dictate the production of doubly charmed 
baryons. We speculate that FSIs play an essential role on the suffering puzzles. Further- 
more, we propose that the FSIs are arisen from the inelastic scatterings D°Df — > 2 C A+ and 
D°D° -> At AT and lead to B -> DDf -> H„At and B -> D°D°K -> AfAzK. The illus- 



trated flavor diagrams are displayed in Fig. ^ For convenience, we will use B 



: A + 

JC c[DD+] 





(a) 



(b) 



FIG. 1: The flavor diagrams for (a) B + — > H[?A+ and (b) B + — > A+A c K + decays, in which the 
former is arisen from B+ D°D+ while the latter is via B + — > D°D°K + . The dotted symbols 
represent the weak interactions. 

and B — > A^A~^ 0l)0 ^K instead of the decaying chains. Since charged B decays have better 
significance, in the following analysis we will concentrate on charged modes. 

To be more clear for the calculations and the effective interactions induced by one-loop 
effects, we replace Fig. ^ of flavor diagrams with Fig. [21 of effective diagrams, where (a) and 
(b) are for two-body and three-body decays, respectively, and the squared (dotted) symbol 
stands for the strong (weak) interactions. Thus, according to Fig. |2fa), the decay amplitude 
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FIG. 2: The effective diagrams for (a) two-body and (b) three-body doubly charmed baryonic 
decays, where double-line denotes the B meson and squared (dotted) symbol stands for the strong 
(weak) interactions. 

for B — > H C A+ could be expressed as 

a^r ■=■ - f A 9DDsB c B c {rn 2 B ) g B DD s (m 2 B ) m 

where u(v) expresses the Dirac spinor field of baryon (antibaryon), e is used for removing 
singularities, mo M is the mass of D( s y c/bdd s (^i?) an d 9dd s b c b c (j^b) denote the effective 
weak coupling for B — D — D s interaction and effective strong coupling for D — D s — B c — B c 
interaction with B c being a charmed baryon, respectively. Since D and D s are the same in 
the SU(3) flavor symmetry except the carrying isospins are different, we will assume that the 
effective strong interactions are the same for D and D s except the effects from the isospin of 
system. Hence, we will use QddBcB c instead of Qdd s b c b c - Similarly, the same assumption is 
also applied to charmed baryons. We note that for two-body decays, because the rest frame 
of charmed baryon-pair is just the rest frame of B meson, we have set the invariant mass 
of charmed baryon-pair, denoted by \/W, to be the mass of B meson. Therefore, the 
values of effective couplings are taken at the m# scale. However, for three-body decays, the 
scale should be chosen at Vk? and is a variable. Although the momentum variable in the 
loop integration in principle has no limit, the main contribution will arise when the on-shell 
condition for the intermediate states D and D s is satisfied, i.e. we can regard the production 
of doubly charmed baryons as the process B — > S C A + = + , as illustrated by Fig.^a). Hence, 

c [UU S J 

in terms of narrow width approximation, l/[(s — M 2 ) 2 + M 2 T 2 M ] m n5(s — M 2 )/MTm, the 
integration of Eq. could be simplified as 

j n n f 1 ~ ggj (2) 

DAPB) J (27T)* q 2 - ml (p B - qf - m 2 Ds 1 ' 

with Pd s being the spatial momentum of D s meson. Since Qbdd s is associated with the 
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decays B — > DDf, we can easily find its relationship to the BRs of B — > DDf by 

BR(B - DDt) = r B ^i\g BDDs (m 2 B )\ 2 . (3) 

By combining Eqs. ((TJ), (J2J) and (jHJ), we find that the BRs for B — > H C A+ decays are related 
to those for B — > DDf by 

BR(B - Z c At) = (l " (mSc ^ mAj2 ) |^ cBc (m|)| 2 5i?(5 - DD+), (4) 

where rriA c (s c ) is the mass of A C (S C ) baryon and px c is the spatial momentum of A c . 

From Eq. P|). we clearly see that BR(B — > 2 C A+) can give a constraint on the parameter 
gDDBcBc{ m B)- However, since the constraint is only suitable at the mg scale, it will not 
help us to understand three-body decays such as B — > A+A~K in which the involving 
invariant mass of the A+A~ system is below the m# scale and the value varies between 2m\ c 
and tub — rax- Moreover, the effects of FSIs should be strongly related to the momenta 
of final state particles, so that the contributions of energetic light baryons are much less 
important than those of slow heavy baryons. In sum, for giving a suitable effective strong 
coupling for gDDB c B c {k 2 ) which could be applied to various values of invariant mass, we 
will model the effective coupling in Lorentz covariant form. In the following, we show our 
way to determine the form of effective coupling. If we regard the scattering DD — > B C B C 
as a t-channel process in which the intermediate particle is a doubly-charmed baryon, the 
scattering amplitude could be approximately described by 

~g 2 (k 2 ) _ , i , „ ~g\k 2 ) _ 

—u Bc r/5 Pd-j 7 7s Pdv Bc ~ u Bc v Bc , 

m D Pd- Pb c - m Bcc m Bcc 

where the fields of D meson have been neglected, g{k 2 ) expresses the coupling D — £> c — B cc , 

m Bcc denotes the mass of exchanged doubly-charmed baryon and the small contributions 

from Pd — p Bc are neglected. Comparing to charmless baryonic decays, we speculate that 

the dominance of FSIs is due to the particles in the doubly charmed baryonic decays carrying 

lower spatial momenta. In order to display the momentum-dependent strong coupling, we 

further parametrize the coupling to be g 2 {k 2 ) = g 2 m Bcc / \\p Bcl — p Bc2 \ A + m 4 ) where p Bcl{2) 

denote the four momenta of charmed baryons, g c is a dimensionless coupling constant and m 

is the effective mass, which is used to remove the singularity when p Bcl = p Bc2 . Accordingly, 

we model the effective coupling for the inelastic scattering DD — > B C B C to be 



Since the mg cc could be estimated by the masses of constituent quarks, the undetermined 
parameters actually are g c and m. 

Similar to the decays B — > H C A^. B£)+ ., we find that the production of two charmed 
baryons in three-body B decays could be arisen from the the inelastic scattering such as 
B — > h+h~, DO p ,K. According to the Fig.^b), the decay amplitude for B — > k+k~K could 
be written as 

J (2vr) 4 q 2 - m 2 D + xe {p B - q) 2 - m 2 D + xe 
= C/ g D DB c Bc(k 2 ) gBKDD(k 2 ,cos9)I D (k 2 )u Ac v Ac , (6) 

where Cj denotes the factor from the isospin wave function of system, taking to be 1/ \pl here, 
and gBKDD{k 2 , cos 9) represents the decay amplitude for B — > D°D°K decays. As known 
that the decay amplitude is Lorentz invariant, for convenience, we will set the working frame 
in the rest frame of charmed baryon-pair. Hence, k 2 stands for the invariant mass of A+A~ 
and is a variable. Angle 9 is the polar angle of A+ with respect to the momentum direction 
of K meson. Since B — > D°D°K are color-allowed processes, it is a good approximation to 
assume that the processes are dominated by factorizable effects. Thus, by employing the 
factorization assumption, the decay amplitude for B — > D°D°K could be written as 

^dd(£ 2 ,cos#) = ^V cb V c * s af(DK\cYs\0) (D\h»c\B) (7) 

where a\ s is the effective Wilson coefficient and the transition matrix elements are defined 

as 



(D(p D )\b lfl c\B(p B )) 
(D(p D )K(p K )\c lfl s\0) 



m B — m n 
[Pb + Pd)v ~ 2 ^ 

m 2 D -m 2 K 
[Pd ~ Pk)^ ^ 



Ff K (Q 2 ^ 
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with q = Pb — Pd and Q = Pd + Pi<- In order to obtain the information of Qbkdd on k 2 
and 9, we have to analyze the decays B — > D°D°K themselves. Since the production of 
doubly charmed mesons is similar to that of doubly charmed baryons, in the following, k 2 
and 9 will be regarded as the invariant mass of D°D° and A+A~ systems and the relative 
angle between Poo t A+ an d Pk, respectively. For deriving the decay rates as a function of 
invariant mass k 2 and angle 9, the coordinates in the k 2 rest frame are chosen as follows: 
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p B {K) = (E B (k),0,0, \p\) with E B(K) = (m 2 B ± k 2 - m 2 K )/2y/W and \p\ = y/E% - m 2 K and 



Pxix) = (v k 2 /2, ±\px\ sin^, 0, ±\pd\ cos9) with \px\ = yk 2 /2- yl — {2mx) 2 /k 2 in which 
X could be _D° meson or A+ baryon. Based on these coordinates, it is obvious that q 2 and 
Q 2 are functions of k 2 and 9. This is the reason why we set the coupling Qbkdd as a function 
of k 2 and cos#. Hence, by including the phase space for three-body decay and using the 
chosen coordinates, the differential decay rate for B — > D°D°K is expressed by 



dT(B — > D°D°K) \p K \ /j2 „. l2 / (2m D ) 2 

with \p K \ = E'k — m 2 K and E' K = {m 2 B — k 2 + m 2 K ) /2m B . In terms of Eq. pj). similarly, 
we can also obtain the differential decay rate for B — > A^A~K decays as 

dT(B -> A+A-AT) fe| ■ . 2 , |2 / _ (2m Ae ) 2 \ 3/2 

dk 2 d cos 9 2 7 7r 3 m|' M Jl V ^ 2 / 

x|^iJD(A; 2 ,cos0)| 2 |C7 JtoJ3Bc B c (A; 2 )| 3 . (10) 

From Eqs. (J7J), © and (fTUj) . we immediately know that if we can determine gBKDD(k 2 , cos 6*) 
from 1? — > D°D°K and gDDB c B c {k 2 ) from I? — > 5 C A+, then in principle we can make predic- 
tions on the BRs of B — > A+A~K. 

In order to understand whether we can predict the BRs for three-body doubly charmed 
baryonic decays, we have to know how many unknown parameters can be controlled by the- 
oretical calculations and parametrizations and constrained by experimental measurements. 
At first, we discuss the unknowns of Eq. (J5J) for Qddb c b c - As mentioned before, m Bcc could 
be estimated by the masses of constitute quarks. If we chose m c = 1.6 GeV and the light 
quark m q = 0.3 GeV, then we get m Bcc = 3.5 GeV. The estimated value is close to the 
particle 5+ (3520), observed recently by SELEX at FNAL Therefore, there remain two 
uncertain parameters g c and rh for effective coupling Qddb c b c - However, the observation 
of B + — > H°A+ only can constrain their ratio. Thus, without further assumption, one un- 
known will be left. Nevertheless, we still can set the value of rh by a proper approach. As 
stated early, the usage of rh is to avoid the suffering problem when the spatial momentum 
of charmed baryon in the k 2 rest frame is approaching vanishment. For two-body decays, 
we can easily know that the spatial momentum of A+ is |p A +| = 1-15 GeV. Therefore, the 
results are insensitive to rh if its value is less than 1.5 GeV. However, for three-body decays, 
the available range of spatial momentum of A+ is |p A +| = [0, 0.71] GeV. Needless to say, the 
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value choice of m will influence on the BRs of B — > A+A~K significantly. To get a balance 
in both kinds of decays, the in can be chosen so that its value will not affect the two-body 
decays and still can protect three-body from overestimation. Since \pb c i ~ Pb c i \ ~ 2|pX c | in 
both two-body and three-body decays, we find the better value for m is around 1.0 GeV 
which is similar to the magnitude of momentum in two-body decays. Hence, if we accept 
this value as input, B + — > 2°A+ can directly give a bound on g c . 

Next, we concentrate on the effective coupling of Eq. ((7|). From the equation, we see that 
the main unknowns are from the form factors i'i(o) and Fff£ . We can use a proper QCD 
approach, such as perturbative QCD [l£|, relativistic quark model and light-front quark 
model [18[ etc, to calculate the transition form factors F^ y According to the results of 
light-front quark model, the explicit expressions can be given by |18] 

. 2 0.67 2 0.67 

l(g ' = (l-1.25g 2 /m| + 0.39(g 2 /m|) 2 )' ° (<? ' = (1 - 0.65q 2 /m%) ' ( ' 

As for the form factors F^, since there exists no good method to estimate the time-like 
form factors, what we can do is as usual to parametrize them to be 



Fi(Q 2 ) = DK % 2 , Fr(Q 2 )= r-^V , (12) 
1 + a u ^Q Ji jm x l + Q z /m x 

where we set m x = ttlb — vtid as the maximum invariant mass of DK system. We note that 
the parametrizations in Eq. (|12j) do not display the behavior predicated by PQCD at the 
large Q 2 , i.e. F DK — > 1/(Q 2 ln(Q 2 /A 2 )) j^]; nevertheless it could be a good approximation 
since the dominant contributions for B — > D°D°K are close to the threshold region of 
the invariant mass in the DK system. In terms of the parametrizations, we see that two 
parameters are needed to be determined. However, so far only the BRs of B — > D°D°K are 
measured. Hence, we only can determine the allowed ratio of a DK and F DK (0) by the data 
of B — > D°D°K. If experiments can provide the information on angular distribution, we 
can further fix the remaining unknown. Nevertheless, we can still obtain some information 
on F DK (0) by considering the two-body decay B° — > D~K + . It is known that the decay 
is g overned by annihilation topology |21( and the BR is measured to be (3.8 ±1.3) x 10~ 5 
|l2j. In terms of factorization assumption, the decay amplitude for B° — > D~K + can be 
expressed by 

m(b°^d;k + ) = ^=v cb v: d aff B (D;K + \s£ B c\o) 

= ^V cb V: d aff B (ml s -m 2 K )F^ K (ml) 



where f B is the decay constant of B meson, a| ff is the color suppressed effective Wilson 
coefficient and nonfactorizable effects have been included, and the time-like form factors 
defined in Eq. (JHJ) are used. If we neglect the differences between D s and D mesons, we 
can get the value of form factor F® K at mg scale. By taking V c b = 0.041, f B = 0.2 
GeV, m Ds(i n = 1.969(0.493) GeV, G F = 1.166 x 10~ 5 GeV~ 2 , t b o = 1.55 x 10~ 12 s Q, 
af = 0.275 |2J, and BR(B° -» D~K + ) = 3.5x 10~ 5 , we get F DX (m|) « F° sK (m 2 B ) = 1.0. 
Relating the extracted value to the Fq K (Q 2 ) parametrized in Eq. ()12j) . we immediately 
obtain F DK (0) = 3.41. If we adopt this value as input, the BRs of B — > D°D°K can fix the 
parameter cr ^. Consequently, by including proper measurements, we can determine the 
introduced parameters with errors from experiments. 

In order to calculate the numerical values, besides the values taken before, we set other 
inputs as follows: m H o = 2.471 GeV, m Ac = 2.285 GeV, m D o = 1.865 GeV, t b + = 1.67 x 
0" 12 s, BR(B + -> D°D+) = (1.3 ± 0.4)%, BR(B+ -> D°D°K) = (1.9 ± 0.4) x 10^ 3 



12( and a\ s = 0.986 [22j. Although in general it is not necessary to regard mg cc = 3.5 
GeV, in = 1.0 GeV and F DK (0) = 3.41 as the inputs, however, it will be interesting 
to see whether by using these values, the predicted BRs of B — > A^A~K are consistent 
with the observations. In addition, we find that actually the BRs of B — > D°D°K are 
insensitive to the value of a DK . That is, in our parametrizations, it is a good approximation 
to take a DK = 1.0. Accordingly, we have BR(B + -> D°D°K + ) = 2.1 x 10~ 3 and the 
differential BR as a function of invariant mass M(D°D°) for B + — > D°D°K + is presented in 
Fig-E^a)- To demonstrate the correlation between g c and BR{B + — > A^A~K + ), we further 
set g c = (3.1, 3.3 , 3.5) so that the BR(B + -> S°A+) = (3.4, 4.4, 5.6) x 10~ 3 and BR(B + -> 
A+A~K + ) = (6.6, 8.5, 10.9) x 10 -4 ; then, we display the corresponding differential BR 
for B + — > A+A~K + as a function of invariant mass M(A+A~) with dash-dotted, solid 
and dashed lines, respectively, in Fig. Efb). From these results, we see clearly that the 
observations of B — > H C A+ and B — > A+A~K could be understood by the similar FSIs. 
Furthermore, to illustrate the influence of the introduced parameters F DK (0), a DK , in and 
g c , we free these parameters and let the data decide the allowed ranges. The results are 
shown in Fig. 0J By Fig. Ufa), we confirm that the BRs of doubly charmed mesonic decays 
are not sensitive to a DK . Taking la error on BK(B + — > D°D°K + ), the range for the value of 
F DK {0) is (3, 3.9). With the bound (2.4, 5.3) x 10" 3 on BR{B + -> E°A+) 0, from Fig.H 
we find that BR(B + — >• S°A+) is sensitive to g c but insensitive to in and we also see that 
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3.6 3.8 4.0 4.2 4.4 4.6 4.8 4.55 4.60 4.65 4.70 4.75 4.80 

(a) M(D°D°)(GeV) (b) M(AV)(GeV) 

FIG. 3: The differential BRs for (a) B + -> D°D°K + and (b) B + -> A+A~K + as a function 
of invariant mass M(D°D°) and M(A+A~), respectively, where we have set F D ^(0) = 3.41, 
cr ^ = 1.0, mB cc = 3.5 GeV and rh = 1.0 GeV. The dash-dotted, sold and dashed lines represent 
the contributions when g c = (3.1, 3.3, 3.5) and their corresponding BRs for B + — > A+A~K + are 
(6.6, 8.5, 10.9) x 10~ 4 , respectively. 




D°D°K + ) and (b) between (g c , rh) and BR{B + -> H°A+). The values (in units of 10" 3 ) in 
the figures (a) and (b) denote the low and upper bounds of current data on BRs. (c) the possible 
BR for B + — ► A^A^Ti^ when the involving parameters are satisfied with the bounds of figures 
(a) and (b), where the values (in units of 10 -4 ) shown in the figure come from the data within la 
statistical and systematic errors. 



a wide range of BR(B + — > A^A~K + ) is allowed, including the data shown. We note that 
the values of BR (in units of 10 -4 ) shown in the figure (c) only contain la statistical and 
systemic errors. By the figure, we also find that the three-body doubly charmed baryonic 
decays are much more sensitive to g c and rh, in particular rh. Clearly, based on our modeling 
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of Eq. (jSJ), the parametrizations of Eq. (|12|) and the limited information on data, in general 
we cannot make concrete predictions on the decays B —> A+A~K. 

Although at present we cannot give precise predictions on doubly charmed baryonic 
decays, however, in terms of the proposed mechanisms we can have some implications on 
other unobserved processes such as B — > H°S C and B — > AfH^K etc. Since and A+ 
have the same quark contents except the former is isospin 1 = 1 while the latter is I = 0, by 
utilizing Eqs. (0) and (0) and introducing a corrected factor Cj = a/2/3 for isospin, we can 
easily obtain the BRs for B — > S°S C . With the values of g c and rh which are satisfied with 
the measurement of B + — > S°A+, the estimated BR(B + — > H°A+) is presented in Fig. Efa) 
in which we have set ms c = 2.452 GeV and the values in the figure denote the available BR 
in units of 10~ 3 . It is clear that although the phase space of the H°£+ mode is less than that 
of the S°A+ mode, by the modeling of Eq. @, we predict that B + — > and B + — > S°A+ 

have similar BRs, even the former could be larger than the latter. Similarly, in terms of 




g s a DK /F DK 

(a) (b) 

FIG. 5: The legends of (a) and (b) are the same as Fig. 0Kb) and (c) but for the decays B + — > 
and B + — > A^T^K + , respectively. 

Eq. (|10J) with the values of the parameters satisfied with the observation of B + — > A+A~K + 
within la error, the estimated BR(B + — > AfYtJK + ) is displayed in Fig. Efb) where the 
values in the figure are the available BR in units of 10 -4 . From the figure, we see that the 
BR of B + — > AfH^K + could be as large as 10~ 4 . If we extend the FSIs to the processes 
associated with the creation of the strange quark-antiquark pair, i.e. ss pair instead of dd 
pair of Fig. ^a), and assume that the relevant strong couplings are the same, interestingly 
we obtain BR(B + — » f2 c H+) pa 5.5%. Although f2 c 5+ has much less allowed phase space, 
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due to the enhancement of Eq. (J5J, we see that the BR for B + — > f2 c H+ is not suppressed. 
However, since we have assumed that the effective couplings of FSIs for the production of 
doubly charmed baryons are the same, one cannot take the result of B + — > f2 c H+ seriously. 
After all, the FSIs for real situations could be different in different processes. What we have 
displayed is that our FSIs could enhance the BR of process which the corresponding phase 
space is suppressed. 

Finally, we make some remarks on other possible FSIs for three-body modes, such as 
B E C A C A+A~K and B -> DD+ A+A~K. As known that the BR(B+ -> A+A~K+) 
is only few factors smaller than BR{B + — > H°A+), if we regard that the FSIs are inelastic 
scattering, by the loop and phase space suppressions, we expect that the contributions 
through B — >■ E C A C — >■ A+A~K should be small. As for the decay B + — >■ D°D+ — >■ A+A~K + 
shown in Fig. Efa), if its contribution is significant, one can speculate that the same FSIs 
could also arise from the doubly charmed mesonic decays such as B + — > D°D°K + , shown 
in Fig. EJb). From Fig. Efa) and (b), we conjecture that the former should be smaller than 
the latter due to one more quark pair being produced. Since we have assumed that the 
dominant effects for B + — > D°D°K + are factorizable contributions, described by Eq. (JZj), 
if the decay chain B + — > D°Df — > A^A~K + is not negligible, in our analysis it could be 
regarded subleading effect. 




(a) (b) 



FIG. 6: Flavor diagrams (a) for B+ -> D°D+ -» A+k~K+ decay and (b) for B + -> D°D+ 
D°D°K + decay. 

In summary, we have studied two-body and three-body doubly charmed baryonic B 
decays. We find that the observed processes could be produced by FSIs B + — > "E e A* 
and B — > A^A~^ 0E)0 ^K. In terms of proper modeling for the effective strong coupling and 
factorization assumption, we find that with the constrained values of the parameters, the 
estimated results could be compatible with the current data. In addition, we also show the 
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implications of FSIs on B + — > S^S^ and B + — > AfHJK + and get that their BRs could be 
as large as BR(B + — > H°A+) and 10 -4 , respectively. In addition, extending the FSIs to 
the processes associated with the creation of ss pair, BR(B + — > Cl c E+) at percent level is 
achievable. 
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